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椭圆 轨道 编队 构 型 的 初始 化 控制 研究 
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1( 中 国 科学 院 空间 科学 与 应 用 研究 中 心 ”复杂 航天 系统 电子 信息 技术 重点 实验 室 ”北京 100190) 
2( 中 国 科学 院 大 学 ”北京 100049) 


摘 要 研究 了 椭圆 参考 轨道 下 的 编队 构 型 初始 化 问题 , 给 出 轨道 面 内 绕 飞 构 型 的 三 种 初始 化 策略 . 基于 状态 转 
移 矩 阵 研究 了 构 型 初始 化 的 两 脉 串 机 动 方 法 ; 基于 高 斯 摄 动 方程 , 以 降低 燃料 消耗 为 原则 , 研究 了 构 型 初始 化 的 三 
迹 向 脉 串 和 四 迹 向 脉冲 方法 , 推导 了 速度 增 量 的 求解 方程 ; 对 上 述 三 种 构 型 初始 化 策略 进行 了 仿真 分 析 , 利用 序列 
二 次 规划 和 遗传 算法 对 计算 结果 进行 优化 . 数值 仿真 表明 , 三 种 初始 化 策略 均 能 完成 构 型 初始 化 , 三 迹 向 脉冲 初始 
化 比 两 脉冲 初始 化 节省 约 41% 的 速度 增 量 , 四 迹 向 脉冲 初始 化 需要 的 速度 增 量 与 三 迹 向 脉冲 初始 化 基本 相等 . 
关键 词 椭圆 参考 轨道 , 构 型 初始 化 , 三 迹 向 脉冲 , 四 迹 向 脉冲 
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雷 博 持 等 : 椭圆 轨道 编队 构 型 的 初始 化 控制 研究 
0 引言 


近年 来 小 卫星 编队 飞行 技术 得 到 快速 发 展 , 与 伟 
统 单 颗 星相 比 , 编队 具备 很 多 优势. 特别 是 随 着 小 卫 
星 、 微 小 卫星 以 及 纳 卫星 技术 的 不 断 进步 , 编队 越 来 
越 受到 重视 , 编队 飞行 技术 已 在 工程 实践 中 得 到 飞行 
验证 [也 

构 型 初始 化 是 航天 器 编队 飞行 中 的 关键 技术 之 
一 已 有 很 多 学 者 对 编队 构 型 初始 化 策略 进行 了 大 
量 研究 , 但 是 大 部 分 研究 集中 于 圆 轨道 编队 初始 化 
问题 , 对 椭圆 轨道 的 初始 化 仍 有 待 进一步 研究 ， 王 
兆 魁 各 研究 了 两 径 向 脉 串 作用 时 , 轨道 面 内 环绕 构 
型 的 初始 化 问题 ; 孟 云 稚 问 在 分 析 单 向 脉冲 对 编队 
的 影响 后 , 研究 了 两 径 向 、 三 迹 向 脉冲 的 初始 化 方法 ; 
张 玉 琨 9 和 朱 户 伟 四 分别 研 究 了 圆 参考 轨道 和 术 
圆 参考 国道 下 编队 构 型 初始 化 的 两 冲 量 问题 

本 文 研究 椭圆 参 考 轨道 下 的 编队 构 型 初始 化 问 
题 . 基于 状态 转移 矩阵 , 对 两 脉冲 初始 化 问题 进行 研 
究 ; 基于 高 斯 摄 动 方程 , 推导 了 三 迹 向 脉冲 和 四 迹 向 
脉 中 的 构 型 初始 化 方法 , 给 出 所 需 脉冲 速度 增 量 的 求 
解 方程 ; 进行 了 仿真 验证 并 对 计算 结果 进行 优化 , 对 
三 种 构 型 初始 化 策略 进行 了 综合 比较 分 析 . 


1 构 型 初始 化 原理 


假设 初始 时 刻 主 从 星 轨道 根 数 近 似 相同 , 在 编队 
构 型 已 知 时 , 根据 构 型 设计 参数 可 以 得 到 从 星相 对 主 
星 的 相对 轨道 根 数 , 一 般 是 轨道 根 数 之 差 . 构 型 初始 
化 是 通过 一 系列 轨道 机 动 控制 , 改变 从 星相 对 主星 的 
轨道 根 数 , 产生 需要 的 相对 轨道 要 素 差 , 进而 形成 所 
需 的 编队 构 型 . 

ur, un, U9 分 别 表示 从 地 心 指向 卫星 方向 、 卫星 
角 动 量 方向 以 及 与 w, un 构成 右手 正 交 坐标 系 的 方 
向 . 设 控制 矢量 = (wr wn uo)' ,wu 对 轨道 要 素 a， 
e i, 0 w, M 的 影响 可 用 高 斯 摄 动 方 程 来 描述 , 有 


da 2a P 
过 本 (esi 十 二 0) 
de 1 

HiFi fu t+lP+r) eosf +relu), 
di 加 7 COSO 

有 
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df rsing 
dt heni™ 
dw 1 


.dn 
a jet-Peos furt (P+r7)sin juo] 一 cos? 
dM 

一 一 一 多 十 [I(Peosf — 2re)ur — (P+7) sin fuel. 
dt he 四 


式 (1) 中 ， 
0=w+f, n=V1l-e?, P=a(l-e’), 
n= Vu/a3, h= VuP. 


从 高 斯 摄 动 方程 可 以 看 出 , 轨道 面 外 的 构 型 初始 
化 相对 简单 . 只 需要 在 相应 位 置 处 施加 法 向 脉冲 , 改 
变 轨道 倾角 或 者 升 交 点 赤 经 , 即 可 完成 轨道 面 外 构 型 
的 初始 化 . 

轨道 面 内 构 型 的 初始 化 比较 复杂 , 需 在 保持 轨道 
半 长 轴 不 变 的 同时 , 改变 相对 偏心 率 、 相 对 近地点 幅 
角 和 相对 平 近 点 角 等 . 本 文 主 要 研究 轨道 面 内 绕 飞 编 
队 构 型 的 初始 化 方法 , 对 经 典 两 脉冲 法 , 重点 对 椭圆 
参考 轨道 编队 初始 化 的 三 迹 向 脉冲 法 和 四 迹 向 脉冲 
法 进行 了 研究 . 


2 构 型 初始 化 策略 


2.1 ”两 脉冲 

根据 两 脉冲 机 动 原理 , 两 脉冲 构 型 初始 化 方法 是 
指 , 给 定 初始 时 刻 to 及 另 一 时 刻 二, 以 及 to 时 刻 从 
星相 对 主星 的 相对 运动 状态 ro, v6，, 在 这 两 个 时 刻 
施加 两 次 脉冲 , 使 去 时 刻 脉冲 作用 后 的 相对 运动 状 
态 为 71, vi, 即 期 望 构 型 所 对 应 的 运动 状态 . 在 该 方 
法 中 , 基于 状态 转移 矩阵 求解 脉冲 速度 增 量 . 其 中 上 
标 --, 二 分 别 表示 脉冲 施加 前 后 . 

根据 椭圆 参考 轨道 相对 运动 动力 学 方程 , 可 以 推 
导出 该 动力 学 方程 的 状态 转移 矩阵 $, 进而 利用 状态 
转移 矩阵 求 出 两 次 脉冲 速度 增 量 的 大 小 站 .其 过 程 
如 下 . 

由 状态 转移 矩阵 可 得 


VD 三 Piz (T1 — p11 * 70), 
v7 = p21 To p22: of (2) 


其 中 状态 转移 矩阵 


p11 p12 
一 ) 3 
021 022 
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进而 可 以 得 到 两 脉冲 矢量 的 大 小 为 


Avo = v0 一 V0， 


Av1 = v1 — V1. (4) 


给 定 一 个 五 ， 根 据 式 (4) 可 以 求 出 一 组 Avo， 
Avi, 对 应 一 个 两 脉冲 机 动 策略 . 另外 , 当 皇 是 未 知 
量 时 , 两 次 脉冲 速度 增 量 之 和 是 时 间 握 的 函数 , 对 这 
一 非 线性 优化 问题 进行 寻 优 , 即 可 得 到 在 给 定时 间 段 
内 两 脉冲 构 型 初始 化 所 需要 的 最 小 速度 增 量 . 

2.2 ”三 迹 向 脉冲 

假设 达到 期 望 构 型 时 , 从 星相 对 主星 的 相对 轨道 
根 数 (轨道 根 数 之 差 ) 为 [Aa, Ae, Ai, AD, Aw,AM]. 
对 于 轨道 面 内 编队 , Ai, A 均 为 零 . 在 理想 情况 下 ， 
为 保证 等 周期 性 旦 不 考虑 半 长 轴 修 正 时 ,Aa 一般 也 
为 零 . 因此 考虑 在 三 个 位 置 处 ， 即 参考 星 真 近 点 角 
为 万 ,， 户 ,， 户 时 分 别 施加 三 次 迹 向 脉冲 以 形成 Ae， 
Auw, AM. 设 这 三 次 迹 向 脉冲 分 别 为 dv1, dv2, dvs. 

从 式 (1) 可 以 看 出 , 利用 迹 向 脉冲 调整 e  w，M 
时 , 均 会 引起 半 长 轴 改 变 . 设 三 次 迹 向 脉冲 所 引起 的 
轨道 根 数 改变 量 为 [ba, 5e, bw, 5M]. 则 根据 式 (1) 可 
以 得 到 


1 
be = 元 tp 十 ri)cos 广 十 7le| duw1 十 


H 72) cos 户 十 roel dv2 十 


(p 
[(p + rs) cos fs + r3e]dvs)}, 
1 


0w = 二 一 [(p 十 71)sin fi dvi 十 
he 


(D 十 r2)sin fo dv (p+ 73)sin fs dvsl], 


0M = 一 (P+) sin 万 dul 十 


(D 十 7a)sin 户 dvz + (p73)sin 户 dvas] 一 
3(Ms jam Mi)ap 3(Ms ER M2 )ap 
en ee i (5) 
式 (5) 中 , 7i;,， Mi 分 别 表示 卫星 在 真 近 点 角 为 fi 处 
的 半径 和 平 近 点 角 . 
对 于 式 (5) 中 的 5A, 第 一 项 根据 高 斯 摄 动 方程 
得 到 , 第 二 项 和 第 三 项 是 脉冲 控制 导致 半 长 轴 发 生变 
化 从 而 引起 的 平均 角 运 动 误差 站 . 第 三 次 脉冲 作用 
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完毕 , 半 长 轴 误 差 再 变 为 零 . 其 推导 过 程 如 下 . 
根据 


n= VS (6) 
可 得 
07 = -了 ja (7) 
又 由 于 
ja = 2 dv, (8) 
则 可 得 
OnN1 一 一 了 du (9) 
同 理 ， 
On2 一 一 了 do (10) 
由 于 
t3—t1= (ms— mi)/n, 
妇 一 妇 = (m3— m2)/n, (11) 
总 的 平均 角 运 动 误 差 为 


Oni(ts x t1) 十 (6n1 十 0n2)(ts AR t2) 一 
Oni(t3 < t1) 十 On2(ts = t2) = 
3(Ms Mi)ap 3(Ms -2 M2 )ap 
加 hri do 和 hro wo (12) 


令 脉冲 引起 的 轨道 根 数 变 化 量 等 于 期 望 构 型 所 
对 应 的 轨道 根 数 差 , 则 可 得 到 如 下 方程 : 
ou = Aa, 
be = Ae, 
0w = Aw, 
SM = AM. (13) 


对 式 (13) 进行 简化 分 析 , 根据 高 斯 摄 动 方程 可 
以 看 出 , 在 分 别 施 加 三 次 迹 向 脉冲 时 , 采用 先 调整 w 
再 调整 e 的 策略 , 在 该 调整 过 程 中 M 同时 得 到 调整 ， 
达到 节省 燃料 的 目的 . 另外 , 为 保证 脉冲 作用 的 效率 
最 高 , 即 ve 前 面 系数 的 绝对 值 最 大 , 可 以 给 定 端点 
值 户 , 户 , 令 户 待定. 

这 样式 (13) 中 包含 4 个 未 知 数 及 4 个 方程 , 可 
以 通过 数值 方法 求解 . 更 进一步 , 当 脉 冲 作 用 点 户 ， 
户 不 固定 时 , 可 以 使 用 序列 二 次 规划 算法 对 及 ，fs 
寻 优 , 先 确定 及, fs, 再 按照 上 述 方法 求解 方程 . 
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2.3 ”四 迹 向 脉冲 
与 2.2 节 类 似 , 考虑 在 4 个 位 置 即 参 考 星 4 个 
不 同 真 近 点 角 对 应 的 位 置 处 分 别 施加 4 次 迹 向 脉冲 ， 
对 e,w, M 进行 调整 , 同时 也 要 保持 Aa 为 零 . 
根据 式 (1) 以 及 2.2 节 内 容 , 略 去 详细 推导 过 程 ， 
可 得 如 下 方程 组 : 
和 


0a 一 dv3 十 


2a2p 
hra 


dva 和 0， 


de = {(p+ri)cosfi+rie] dul 十 


Dp 十 72) cos 户 十 7r2el dv2 十 


( 
[(p 要 国 73) COS fs 十 73e| dvs 十 
( 


[ pT r4) COS fa 十 7T46 dva} 一 Ae, 


= 一 [二 7m)sn 户 dv 村 


(D 十 7r2)sin 户 dvu2z 十 


十 rs)sin f3 dvs+ 


(2 十 7r4)sin fa dvu4] = Aw, 


0M = 一 [+ ri) sin fidvi+t 


(p 十 72)sin 户 du2 十 


(十 rs) sin fadvs 二 (2 十 ra)sin fadva|— 


3(M4 ss MD 2 3(Ma ht Me 
hri hro 

se 一 AM (14) 
hrs 


对 于 式 (14), 当 四 次 迹 向 脉冲 的 作用 点 都 已 知 , 即 广 ， 
户 , 户 , fa 已 知 时 , 该 方程 是 个 四 元 线性 方程 组 , 若 系 
数 和 矩阵 可 道 , 则 可 以 求解 得 到 dv1，dv2, dvs, dva 的 
解析 解 . 

下 面 分 析 脉 冲 作用 点 对 应 的 参考 星 真 近 点 角 的 
选取 . 仍 采用 先 调整 w 再 调整 e, 间接 调整 M 的 策 
略 , 以 节省 燃料 , 同时 也 要 保证 单 次 脉冲 的 作用 效率 
最 高 . 综 上 , 可 以 先 选取 一 组 真 近 点 角 为 [90 270 360 
180], 解 出 一 组 可 行 解 . 

进一步 , 利用 遗传 算法 对 4 个 脉冲 作用 点 进行 
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优化 , 在 优化 的 基础 上 求解 式 (14), 将 优化 结果 与 上 
述 可 行 解 进行 比 对 . 


3 仿真 分 析 


根据 编队 构 型 设计 方法 , 针对 任务 需求 , 可 以 得 
到 轨道 面 内 绕 飞 队 形 在 初始 时 刻 的 主星 轨道 根 数 以 
及 期 望 构 型 所 对 应 的 相对 轨道 根 数 , 结果 列 于 表 1. 

下 面 针 对 第 2 节 的 理论 公式 推导 进行 仿真 验证 ， 
先 在 MATLAB 中 进行 数值 计算 , 解 算出 脉冲 速度 增 
量 , 再 在 STK 中 进行 可 视 化 仿真 分 析 . 

3.1 ”两 脉冲 

基于 2.1 节 的 分 析 可 以 看 出 , 对 于 两 脉冲 构 型 初 
始 化 策略 , 给 定 两 个 时 间 点 , 即 可 根据 解析 解 的 表达 
式 直接 求 出 两 次 脉冲 速度 增 量 大 小 .这 里 采用 对 脉 
冲 作用 点 寻 优 的 方法 来 研究 这 一 问题 , 利用 序列 二 次 
规划 四 对 这 一 非 线性 优化 问题 进行 求解 . 

在 一 个 轨道 周期 内 , 假定 初次 脉冲 作用 点 对 应 时 
刻 to = 0, 仅 对 第 二 次 脉冲 作用 的 时 间 所 进行 寻 
优 ( 单 点 优化 )， 以 此 得 到 第 一 次 脉冲 作用 时 间 点 固 
定时 , 两 脉冲 构 型 初始 化 所 需要 的 最 小 速度 增 量 , 结 
果 列 于 表 2. 

表 2 中 速度 增 量 矢量 在 主星 LVLH 坐标 系 下 (下 
文 同 ). 主星 LVLH 坐标 系 即 主星 质心 轨道 坐标 系 , 坐 
标 系 原点 在 飞行 器 质心 , XY 轴 在 卫星 轨道 平面 内 ， 
X 轴 从 地 心 指向 飞行 右 质 心 方向 , Y 轴 垂 直 X 轴 指 
向 航 迹 方向 , 2 轴 垂 直 于 卫星 轨道 平面 并 满足 右手 螺 
旋 定 则 . 

从 表 2 可 以 看 出 完成 构 型 初始 化 所 需 总 的 脉冲 


表 1 主星 轨道 根 数 与 相对 轨道 根 数 
‘Table 1 Orbital elements of main satellite and 


the relative orbital elements 


轨道 根 数 主星 相对 轨道 根 数 
a/km 7378.137 0 

e 0.0678 一 0.000 336 
i/(°) 90 0 

0Q/(°) 112 0 

w/(°) 0 0.498 439 
MJ/(°) 0 一 0.498 439 
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速度 增 量 大 小 为 4568m's 上 且 第 二 次 脉冲 作用 时 
间 妇 为 3118.542s. 

在 一 个 轨道 周期 内 , 同时 对 两 个 脉冲 作用 点 对 应 
的 时 刻 to, 五 进行 优化 (两 点 优化 ), 结果 列 于 表 3. 从 
表 3 可 以 看 出 , 两 变量 同时 优化 时 , 完成 初始 化 所 需 
速度 增 量 大 小 为 4295m's 两 次 脉冲 作用 时 间 分 
别 为 860.621s 和 5184.706 8. 

对 比 表 2 与 表 3. 与 单 点 优化 相 比 , 两 点 同时 优 
化 相当 于 主星 与 从 星 一 起 运行 到 如 时 刻 , 再 进行 第 
一 次 脉冲 机 动 , 开始 分 离 . 两 个 时 间 变 量 同 时 优化 时 ， 
构 型 初始 化 消耗 的 速度 增 量 小 于 单个 时 间 变 量 优化 
的 结果 . 

利用 优化 所 得 数据 , 在 STK 中 采用 Propagator 
模块 进行 仿真 分 析 ， 这 里 仅 对 两 点 优化 结果 进行 仿 
真 分 析 , 结果 如 图 1~4 所 示 . 其 中 , 图 1 和 图 2 分 
别 表示 单 变量 优化 和 双 变 量 优 化 的 构 型 初始 化 过 程 ， 
图 3 和 图 4 分 别 为 上 述 两 种 情况 下 主 从 星 的 相对 距 
离 在 半天 之 内 的 变化 情况 . 


图 1 单 变量 优化 两 脉冲 构 型 初始 化 过 程 


Fig.1 Formation capture based on two impulses 


maneuver using single-point optimization 


表 2 
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从 图 1~4 可 知 , 两 脉冲 方法 能 够 实现 椭圆 参考 
轨道 下 的 轨道 面 内 构 型 初始 化 ， 而 且 构 型 能 够 稳定 
环绕 . 


图 2 两 变量 优化 两 脉冲 构 型 初始 化 过 程 


Formation capture based on two impulses 


Fig. 2 


maneuver using two-point optimization 


relative distance / km 


time/h 


图 3 单 变量 优化 时 半天 内 主 从 星相 对 距离 变化 


Evolvement of the relative distance between the 


Fig.3 
main satellite and fellow satellite in half a day for 


single-point optimization 


两 脉冲 单 点 优化 的 结果 


Table 2 Results of single-point optimization for two impulse control 


脉冲 作用 时 刻 t/s 速度 增 量 /(m.s ) 
0 [—0.7247 -0.6114 0] 
3118.542 [3.5543 0.6877 0] 


总 速度 增 量 /(m-s-!) 
4.568 


初始 化 所 需 时 间 /s 
3118.542 


表 3 ”两 脉冲 两 点 优化 的 结果 


Table 3 Results of two-point optimization for two impulse control 


脉冲 作用 时 刻 t/s 。 ”速度 增 量 /(m.s”) 
860.621 [—3.6707 0.7368 0] 
5184.706 [~0.00392 一 0.5496 0] 


总 速度 增 量 /(m.s“) ”初始 化 所 需 时 间 /s 


5184.706 


雷 博 持 等 : 椭圆 轨道 编队 构 型 的 初始 化 控制 研究 


3.2 ”三 迹 向 脉冲 
当 式 (14) 中 对 应 的 及 = 7/2, fs = 37, 所 待 
定时 , 可 以 用 差分 的 方法 迭代 求解 该 方程 组 ， 结 果 
列 于 表 4， 从 表 4 可 以 明显 看 出 , 第 二 次 脉冲 作用 
点 的 真 近 点 角 为 289.056°, 需要 的 总 脉冲 速度 增 量 
为 3.012 m.s-1. 

利用 表 4 中 的 数据 进行 仿真 分 析 , 结果 如 图 5 和 
图 6 所 示 . 从 图 5 和 图 6 可 以 看 出 , 施加 三 次 迹 向 脉 
冲 的 方法 能 够 实现 构 型 初始 化 , 而 且 所 得 到 的 绕 飞 构 
型 能 够 长 期 稳定 . 

进一步 利用 序列 二 次 规划 算法 对 三 脉冲 初始 化 
进行 优化 ， 先 对 脉冲 作用 的 两 个 端点 及 ，f3 进行 
优化 ,， 再 在 此 基础 上 求解 方程 . 优化 计算 结果 表明 ， 
当 及 二 27/3, f3 = 8r/3 时 , 所 求解 出 的 三 迹 向 脉 
冲 消耗 速度 增 量 较 小 , 为 2.501 ms-!, 如 表 5 所 示 . 

对 比 表 4 与 表 5 可 知 , 经 过 端点 优化 后 , 完成 三 
迹 向 构 型 初始 化 所 消耗 的 速度 增 量 相 比 端点 固定 时 
可 节省 约 17%; 对 比 表 3 与 表 5 可 知 , 三 迹 向 构 型 初 
始 化 策略 相 比 两 脉冲 构 型 初始 化 策略 可 节省 约 41 兄 
的 速度 增 量 . 


relative distance /km 
a 


time/h 


图 4 两 变量 优化 时 半天 内 主 从 星相 对 距离 变化 


Fig.4 Evolvement of the relative distance between the 


main satellite and fellow satellite in half a day for 


two-point optimization 


表 4 
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3.3 ”四 迹 向 脉冲 

当 四 迹 向 脉冲 的 作用 点 所 对 应 主星 真 近 点 角 依 
次 选取 为 90 270 360 180] 时 , 式 (14) 是 一 个 四 元 
线性 方程 组 , 求解 该 方程 , 可 得 四 迹 向 脉冲 构 型 初始 
化 策略 所 需 脉冲 速度 增 量 大 小 , 结果 列 于 表 6. 

从 表 6 可 以 看 出 , 在 4 个 位 置 处 施加 4 次 迹 
向 脉冲 进行 构 型 初始 化 ,所 需要 的 总 脉冲 速度 增 量 
为 3.399 ms-1. 


图 5 


Formation capture based on three in-track 


三 迹 向 脉冲 构 型 初始 化 
Fig.5 


impulses maneuver 


一 
[4 


oo 


relative distance /km 


全 个 


time /h 


图 6 半天 内 主 从 星相 对 距离 变化 


Evolvement of the relative distance between the 


Fig.6 


main satellite and fellow satellite in half a day 


f1 = 二 7/2, fs 二 37 时 求解 结果 


Table 4 Results when fi = 7/2, fs = 37 


脉冲 作用 点 f/(*) ”速度 增 量 /(m.s”) 
90 [0 0.803 0] 
289.056 [0 —1.465 0] 
540 [0 0.744 ol 


总 速度 增 量 /(ms 一 ) 
3.012 


初始 化 所 需 时 间 /s 
9469.42 
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利用 所 计算 的 数据 进行 构 型 初始 化 仿真 分 析 , 结 
果 如 图 7 和 图 8 所 示 . 从 图 7 和 图 8 可 以 看 出 , 通过 
在 迹 向 施加 四 次 脉冲 的 方法 , 可 以 实现 椭圆 参考 轨道 
下 的 轨道 面 内 构 型 初始 化 , 而 且 能 够 稳定 环绕 . 

对 四 迹 向 脉冲 构 型 初始 化 的 4 个 脉冲 作用 时 刻 
进行 优化 , 利用 遗传 算法 求解 这 一 问题 , 所 得 结果 见 
图 9 和 表 7. 


图 7 四 迹 向 脉冲 构 型 初始 化 


Fig.7 Formation capture based on three in-track 


impulses maneuver 
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从 图 9 和 表 7 可 以 看 出 ,对 脉冲 作用 点 进行 
优化 后 ， 四 迹 向 脉冲 构 型 初始 化 所 需 速度 增 其 大 小 
为 2.499m.s““， 仿真 结果 表明 , 该 数据 可 以 实现 构 
型 初始 化 , 而 且 能 够 稳定 环绕 . 
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图 8 半天 内 主 从 星相 对 距离 变化 


Fig.8 Evolvement of the relative distance between the 


main satellite and fellow satellite in half a day 


表 5 优化 结果 


Table 5 Optimization results 
脉冲 作用 点 f/(*) ”速度 增 量 /(ms~") ”总 速度 增 量 /(m:s”) ”初始 化 所 需 时 间 /s 
120 [0 0.607 0] 2.501 8295.4 
299.586 [0 一 1.208 0] 
480 [0 0.685 0] 
表 6 ”四 迹 向 脉冲 求解 结果 
Table 6 Solutions of four in-track impulses 
脉冲 作用 点 f/(*) ”速度 增 量 /(m-s“") ”总 速度 增 量 /(m:s~") 初始 化 所 需 时 间 /s 
90 [0 0.897 ol 3.399 9469.41 
270 [0 一 1.275 0] 
360 [0 一 0.382 0] 
540 [0 0.844 0] 
表 7 四 迹 向 脉冲 优化 求解 结果 
Table 7 Optimization solutions of four in-track impulses 
脉冲 作用 点 f/(*) ”速度 增 量 /(ms~") ”总 速度 增 量 /(m:s”) ”初始 化 所 需 时 间 /s 
297.344 [0 —0.621 0] 2.499 11619.9 
476.408 [0 0.110 ol 
483.161 [0 1.182 0] 
656.066 [0 一 0.586 0] 
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3.4 燃 耗 对 比 

表 8 列 出 了 以 上 三 种 构 型 初始 化 策略 需要 消耗 
的 速度 增 量 以 及 所 需 时 间 的 综合 比较 . 从 表 8 可 以 得 
出 如 下 结论 . 

(1) 三 种 构 型 初始 化 策略 均 能 完成 顶 圆 轨道 平面 
内 编队 的 初始 化 , 其 中 两 脉冲 方法 消耗 的 速度 增 量 最 
多 , 四 迹 向 脉冲 方法 消耗 的 速度 增 量 最 少 . 

(2) 三 迹 向 脉冲 初始 化 策略 相 比 两 脉冲 初始 化 策 
略 可 节省 约 41% 的 速度 增 量 , 四 迹 向 脉冲 初始 化 策 
略 需要 的 速度 增 量 与 三 迹 向 脉冲 初始 化 策略 基本 相 
等 . 

(3) 本 文 给 出 的 优化 方法 是 有 效 的 . 对 于 两 脉冲 
初始 化 , 优化 后 速度 增 量 节省 了 约 6.7%; 对 于 三 迹 向 
脉冲 初始 化 , 优化 后 速度 增 量 节省 了 约 17%; 对 于 四 
迹 向 脉冲 初始 化 , 优化 后 速度 增 量 节省 了 约 26.5%. 

ee 
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图 9 四 切 向 脉冲 遗传 算法 优化 结果 
Fig.9 ”Optimization results of four in-track impulses 


Using genetic algorithm 


表 8 三 种 初始 化 策略 所 消耗 的 速度 增 量 及 时 间 
Table 8 ”Velocity increment and time of the 


three approaches for formation initialization 


构 型 初始 化 策略 总 速度 增 量 /ms ”所 需 时 间 /s 
两 脉冲 ( 单 变量 优化 ) 4.568 3118.54 
两 脉冲 ( 双 变 量 优 化 ) 4.259 5184.71 
三 迹 向 脉冲 (端点 固定 ) 3.012 9469.42 
三 迹 向 脉冲 (端点 优化 ) 2.501 8295.4 
四 迹 向 脉冲 (无 优化 ) 3.399 9469.41 
四 迹 向 脉冲 (优化 ) 2.499 11 619.9 
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研究 了 椭圆 参考 轨道 下 轨道 面 内 绕 飞 编队 构 型 


的 初始 化 问题 , 推导 了 两 脉冲 、 三 迹 向 脉冲 和 四 迹 向 
脉冲 的 初始 化 方法 , 进行 了 数值 仿真 , 并 对 仿真 结果 
进行 优化 计算 分 析 . 结果 表明 , 两 脉冲 、 三 迹 向 脉冲 
和 四 迹 向 脉冲 方法 均 可 以 实现 椭圆 参考 轨道 平面 内 
编队 构 型 的 初始 化 ， 对 这 三 种 构 型 初始 化 方法 燃料 
消耗 的 对 比分 析 表 明 , 三 迹 向 脉冲 初始 化 策略 相 比 两 
脉冲 初始 化 策略 节省 约 41% 的 速度 增 量 , 四 迹 向 脉 
冲 初始 化 策略 需要 的 速度 增 其 与 三 迹 向 脉冲 初始 化 
策略 基本 相等 . 本 文 研究 结果 可 为 椭圆 轨道 编队 飞行 
任务 的 初始 化 提供 依据 和 参考 . 
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